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First results on charm quarkonia production in heavy ion collisions at the Relativistic Heavy 
Ion Collider (RHIC) are presented. The yield of J/ip's measured in the PHENIX experiment via 
electron-positron decay pairs at mid-rapidity for Au-Au reactions at ysjvjv = 200 GeV are analyzed 
as a function of collision centrality. For this analysis we have studied 49.3 million minimum bias Au- 
Au reactions. We present the J/tp invariant yield dN/dy for peripheral and mid-central reactions. 
For the most central collisions where we observe no signal above background, we quote 90% confi- 
dence level upper limits. We compare these results with our J/ip measurement from proton-proton 
reactions at the same energy. We find that our measurements are not consistent with models that 
predict strong enhancement relative to binary collision scaling. 
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I. INTRODUCTION 

Lattice Quantum Chromodynamics (QCD) calcula- 
tions indicate that there is a transition of nuclear matter 
from confined to deconfined quarks and gluons at a tem- 
perature of order T c = 170 MeV. Characteristic of this 
deconfined state of matter is the dynamic screening of the 
long-range confining potential of QCD. Color screening is 
predicted to reduce the attraction between heavy quark- 
antiquark pairs, and thus leads to a decrease in the ratio 
of hidden charm and beauty (quarkonia) to open charm 
and beauty Q, Thus, one expects a suppression of 
quarkonium states ( J/ip, ip', Xc T(ls, 2s, 3s)) depending 
on their binding energy and the temperature of the sur- 
rounding system. 

In relativistic heavy ion collisions a state of a de- 
confined thermalized quark-gluon plasma may be cre- 
ated. Measurements in Pb-Pb reactions at t/snn = 17.3 
GeV by the NA50 experiment Q show a suppression of 
heavy quarkonia production relative to "normal" nuclear 
absorption, the dissociation of cc pairs by interactions 
with the nucleons into separate quarks that eventually 
hadronize into D mesons 0, ID, HI This suppression 
has been interpreted in the context of color screening in a 
quark-gluon plasma El Q , additional absorption with co- 
moving hadrons 0, E| i an d multiple scattering between 
the charm quarks and the surrounding medium |lll IT^ . 

At RHIC energies, where of order 10 cc pairs are pro- 
duced in central Au- Au reactions 0, 0] , some models 
predict an enhancement of heavy quarkonia due to cc co- 
alescence in a quark-gluon plasma |l5j . detailed balance 
of D + D <-> J/ip + X 0, and/or statistical J/ip pro- 
duction ^tJ- I n addition, at RHIC energies initial state 
effects of shadowing and possible parton saturation may 
play a role in initial charm production [Tsj . Disentangling 
these competing effects will require a systematic study 
of yields of various quarkonium states in different col- 
liding systems (proton-proton, proton (or deuteron)-ion, 
and ion-ion) and over a wide kinematic range in terms of 
transverse momentum and xf- 

We report here the first results on J/ip production via 
electron-positron decay pairs at mid-rapidity from Au- 
Au collisions at v / saT/v = 200 GeV from data taken dur- 
ing Run- 2 at RHIC in 2001. For peripheral and mid- 
central Au-Au collisions, we present the most probable 
yield values, while for central reactions, we observe no 
signal above background and thus quote 90% confidence 
level upper limits on J/ip production. 



II. PHENIX EXPERIMENT 

The PHENIX experiment is specifically designed to 
make use of high luminosity ion-ion, proton-ion, and 
proton-proton collisions at the Relativistic Heavy Ion 
Collider to sample rare physics probes including the J/ip 



and other heavy quarkonium states. The PHENIX exper- 
iment includes two central rapidity spectrometer arms, 
each covering the pseudo-rapidity range \rj\ < 0.35 and 
an interval of 90 degrees in azimuthal angle (p. The 
spectrometers are comprised from the inner radius out- 
ward of a Multiplicity and Vertex Detector (MVD), Drift 
Chambers (DC), Pixel Pad Chambers (PC), Ring Imag- 
ing Cerenkov Counters (RICH), Time-of-Flight Scintil- 
lator Wall (TOF), Time Expansion Chambers (TEC), 
and two types of Electromagnetic Calorimeters (EMC). 
This combination of detectors allows for the clean iden- 
tification of electrons over a broad range in transverse 
momentum. Further details of the detector design and 
performance are given in |19|. 

The Au-Au event centrality is estimated using the com- 
bined data from our Beam-Beam Counters (BBC) and 
Zero Degree Calorimeters (ZDC). While the ZDCs mea- 
sure forward neutrons that result from fragmentation of 
the colliding nuclei, the BBCs are sensitive to charged 
particles produced in the collisions. Together, both de- 
tectors yield information on the impact parameter of the 
nuclear reaction |20|. These observables, combined with a 
Glauber model for the nuclear geometry, allow us to de- 
termine different collision geometry categories, referred 
to as centrality ranges [2l|. 

For the analysis presented here, the electron (positron) 
momentum and charge sign is determined from track- 
ing using the DC and the PC and then projecting back 
through the PHENIX axial magnetic field to the colli- 
sion point determined by the BBC 22] . The momen- 
tum resolution achieved is Sp/p = 0.7% © 1.0% x p (in 
GeV/c). Electrons are cleanly separated from the large 
background of charged pions and kaons by associating the 
tracks with at least three active photomultiplier tubes in 
the RICH 23] . In addition, we compare the track mo- 
mentum (p) to the energy (E) measured in the Electro- 
magnetic Calorimeter. The E/p ratio is used to further 
reduce the pion contamination in the electron sample. 
Pions typically deposit only a fraction of their energy in 
the calorimeter whereas electrons deposit all of their en- 
ergy. These selections are augmented by requiring that 
the calorimeter shower position and time-of-flight infor- 
mation agree with the track projection. Thus, we obtain 
a clean sample of electron and positron candidates with 
less than 5% contamination. 



III. DATA SELECTION AND TRIGGERS 

The Au-Au data at ^/sJTn — 200 GeV used in this 
analysis were recorded during Run 2 at RHIC in the fall 
of 2001. For our "minimum bias" Au-Au event selec- 
tion, we use a Level- 1 trigger that requires a coincidence 
between our BBCs. We place an additional offline re- 
quirement of at least one forward neutron in each of our 
ZDCs to remove beam related backgrounds. Our "min- 
imum bias" sample includes 92% of the 6.9 barn Au-Au 
inelastic cross section [2l| . We further restrict our analy- 



4 



sis to 90% of the inelastic cross section to remove a small 
remaining contribution from beam related background 
events. 

We observed a Au-Au inelastic collision rate that in- 
creased during the running period from 100 to 1200 
Hz. The Level-2 triggers are implemented in a personal 
computer-based farm with 30 processors in Run 2, as part 
of the PHENIX Event Builder The Level-2 J/ip trig- 
ger algorithm identified electron candidates by starting 
with rings in the RICH and then searching for possible 
matching showers in the EMC. The EMC search window 
based on the RICH ring is obtained from a lookup ta- 
ble, generated using Monte Carlo simulations of single 
electrons. Possible matches were assumed to be electron 
candidates, and the electron momentum was taken to be 
the EMC shower energy. The invariant mass was cal- 
culated for all electron candidate pairs within an event, 
regardless of the candidate's charge sign. If the invariant 
mass was higher than 2.2 GeV/c 2 , the pair was accepted 
as a J/ip candidate, and the entire event was archived. 
The Level-2 trigger provided a rejection factor of order 30 
relative to our "minimum bias" Level-1 trigger sample. 

An additional offline requirement was imposed that the 
collisions have a z- vertex satisfying \z\ < 30 cm in or- 
der to eliminate collisions taking place near the PHENIX 
magnet. After this selection, we have analyzed 25.9 mil- 
lion "minimum bias" Au-Au reactions as triggered by our 
BBC Level-1 Trigger. In addition, from the high lumi- 
nosity period of running, we also processed 23.4 million 
"minimum bias" events with our J/ip Level-2 trigger. 



TABLE I: Statistical results for J/ip counts are shown for 
three exclusive centrality ranges. Shown are the number of 
unlike and like sign counts within the mass window (2.8 < 
m < 3.4 GeV/c 2 ). Also shown are the most likely signal value 
with the 68% statistics confidence interval (for the peripheral 
and mid-central cases), and the 90% confidence level upper 
limits. 
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Unlike Sign 
Counts 


Like Sign 
Counts 


Most Likely 
Signal 


90% 

C.L. 


00-20% 


33 


41 





9.9 


20-40% 


16 


8 


n+4.8 
8 -4.1 


14.4 


40-90% 


7 


2 


c+3.1 
J -2.6 


9.3 



formation to be utilized in the exact distribution of the 
candidates within the mass window. However, we have 
found that this does not add to the significance of the 
result given the low counts and the lack of constraint on 
the J/ip and background line shape. 

Table^shows the number of unlike and like sign counts 
within the mass window. For given observed counts of iV; 
(like sign) and N u (unlike sign), the likelihood L(vi,v u ) 
for the expectation values V\ and v u is given as 



L(v u i> u ) 



Ni - 



NJ 



(1) 



We then integrate L{i>i,v u ") to give the likelihood L(v s ) 
for the expectation value of the net signal counts v s — 

v u - n 



IV. J/ip SIGNAL COUNTING 



L{v s ) 



L(v u v u )8(v s -i> u + i>i)dvidv u (2) 



For three exclusive centrality bins, 0-20%, 20-40%, and 
40-90% of the total Au-Au cross section, we show the di- 
electron invariant mass distributions for unlike sign pairs 
(e + e~), like sign pairs (e + e + or e~e~) and the subtracted 
difference in Figure^ The number of J/ip counts for each 
centrality range is determined from the number of signal 
counts above "background" within a fixed invariant mass 
window. The PHENIX acceptance and Level-2 trigger ef- 
ficiencies are the same within a few percent for unlike sign 
pairs and like sign pairs in the J/ip mass region. There- 
fore, the sample of like sign pairs is a good representation 
with no additional scale factor of the "background" due 
to simple combinatorics. 

In order to extract a J/ip signal strength, we employ a 
counting method where we subtract from the number of 
unlike sign pairs the number of like sign pairs in the mass 
window 2.8 < m < 3.4 GeV/c 2 . We have chosen a wide 
invariant mass window to be consistent with the si gna l 
extraction method from our proton-proton analysis |24| . 
and to limit our sensitivity to the exact mass width value. 
Although we expect from our Monte Carlo studies a mass 
width of order 60 MeV, we cannot quantitatively verify 
this even with our proton-proton data sample due to low 
statistics. We note that in principle there is more in- 



We show the likelihood distribution L(y s ) for the mid- 
central (20%-40%) events in Figure^ Although L(u s ) is 
normalized such that L(y s )dv s = 1, it has a non-zero 
probability for negative expected net signal value (y s < 
0). Since the unlike sign contains signal + background, 
and the like sign contains only background, the only 
physically allowed values for v s are greater than or equal 
to zero. Thus, we remove the probability range corre- 
sponding to v s < 0, and re-normalize the remaining prob- 
ability integral to one [25| , as shown in the same Figure|21 
We then determine for each centrality the 90% confidence 
level upper limit, and the 68% confidence interval around 
the most likely value for the peripheral and mid-central 
ranges. These values are shown in Table [I] 

Since the net signal is negative for the 0-20% central 
event class, we can only quote a 90% confidence level 
upper limit. Also, even for the 20-40% and 40-90% cen- 
trality classes, the signal observed is not significant at 
the two standard deviation level and thus we also show 
90% confidence level upper limits for completeness. The 
limited statistical significance of the results is clear from 
the mass distributions shown in Figure ^ 

In the intermediate mass region below the J/ip, 2.0 < 
m < 2.8 GeV/c 2 , the shapes and absolute yield of like 
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FIG. 1: (Color online) Dielectron invariant mass distribution in Au-Au reactions (top row: most central, 0-20% central, 
middle row: mid-central, 20-40% central, and bottom row: peripheral, 40-90% central) for unlike sign pairs containing sig- 
nal+background (left column), like sign pairs containing only background (center column) and the subtracted difference (right 
column) . 
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FIG. 2: (Color online) The Poisson statistical likelihood 
distribution as a function of the expected net signal. The 
distributions are for the mid-central case of N U nlike=l6 and 
Nuk e —8. The dashed curve is the likelihood distribution, and 
the black is after eliminating the unphysical net signal less 
than zero and re-normalizing. Vertical lines are shown to in- 
dicate the most likely value (8), the 68% confidence interval 
values, and the 90% confidence level upper limit. 



contribution is less than the statistical uncertainty of the 
data (2|| • A complete analysis of the intermediate mass 
region will be presented elsewhere [22} . 

Our method of measurement does not separate contri- 
butions from decay feed-down from other states such as 
the Xc and ip' , and thus our J/ip counts include these feed 
down sources. Our acceptance and efficiency is identical 
for the resulting J/ip from these decays as for prompt 
J/ip, and thus they enter our signal weighted simply 
by their relative production and branching fraction into 
J/ip + X. Another contribution may result from B me- 
son decays into J 'ftp. However, if we assume a bb cross 
section in proton-proton reactions a b ^ sa 2 — 5 fib psj 
and that beauty production scales with binary collisions, 
we would expect a contribution of order 1-4% relative to 
primary J/ip or J ftp from \ c . This percentage contribu- 
tion calculation assumes that the primary J/ip also scales 
with binary collisions. If primary J/ip are substantially 
suppressed, the B meson decay contribution would con- 
stitute a larger fraction of our measured J/ip, especially 
at higher px ■ 

It should also be noted that some signal in the un- 
like sign pairs from correlated open charm cc — > D(—> 
e + + X) + D(—> e~ + X) will contribute in our mass 
window. Assuming binary scaling in charm production 
with a proton-proton cross section a C c ~ 650 /ib fLU ] , this 
contribution in the J/ip mass region is estimated to be 
about 0.1 events in the 0-20%, 0.05 events in the 20-40%, 
and 0.02 events in the 40-90% centrality bins. 



sign and unlike sign dielectron mass distributions are well 
reproduced by an event mixing method within a few per- 
cent. This indicates that most of the dielectron pairs 
are from uncorrelated electron and positron candidates. 
They are originating from Dalitz decays, photon conver- 
sions, open charm/beauty semi-leptonic decays, and a 
small contamination of mis-identified hadrons. The un- 
like sign pairs should also have a component from semi- 
leptonic decays of charm and anti-charm pairs, but its 



V. J/ip YIELD CALCULATION 

We quote our results as the branching fraction of 
J/ip -> e+e" (B=5.93±0.10x 10" 2 [H) times the invari- 
ant yield at mid-rapidity dN / 'dy\ y —o. We calculate this 
quantity for three exclusive centrality ranges as detailed 
below. 



B^ = ^ -x-Lx I (3) 

dy N m b- e vt + (Qvl2-eff X Nl v l2-evt) ^acc-eff X e cent 

I 



The number of signal counts Nji^ from both the 
"minimum bias" and Level-2 triggered event samples are 
shown in Table [I] The number of events from the "mini- 
mum bias" sample N m i,_ evt is 25.9 million Au-Au events. 
The number of effective events sampled by the Level- 
2 trigger is (e;„/ 2 - e // X iVjuja-eut)) which is the Level-2 
trigger efficiency times the number of events processed 
by the Level-2 trigger, 23.4 million Au-Au events. This 
formulation appropriately weights the two data samples 



by the expected number of J/ip. 

The efficiency of the Level-2 trigger £i v i2-eff was de- 
termined by running the trigger algorithm on simulated 
J/ip and carrying out a full offline reconstruction of the 
resulting electron-positron decay pair. The efficiency was 
calculated via counting the fraction of successfully recon- 
structed J/ip events that were also found by the trigger. 
In these trigger simulations, the channel-by-channel cal- 
ibrations for the RICH and EMC were used to convert 
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the simulated signals into realistic values representative 
of a specific period in the run, before passing them to the 
Level-2 trigger. 

The overall J/tp trigger efficiency from the trigger sim- 
ulations was €[ v i2- e ff — 0.75 ± 0.04. The systematic er- 
ror was determined by studying the dependence of the 
trigger efficiency on the collision vertex position, the as- 
sumed J/tp transverse momentum and rapidity distribu- 
tion, collision centrality, and the period of the run from 
which the channel-by-channel calibrations were taken. 
After evaluating all of the above dependencies, we assign 
a 5% systematic error to the J/tp trigger efficiency. 

The efficiency result from the trigger simulations was 
confirmed using real data in two ways. First, the mini- 
mum bias data sample in Au-Au collisions was analysed 
to calculate the J/tp trigger efficiency. This was done by 
taking the ratio of the events that have an electron pair 
with invariant mass between 2.8 and 3.4 GeV/c 2 that 
fired the J/tp trigger to all of the events having electron 
pairs in that invariant mass range. The trigger efficiency 
estimate from this check is 0.67 ± 0.10 (stat). Second, 
the triggers were run on a sample of 26 events from the 
proton-proton data set that passed all of the J/tp cuts in 
the offline analysis and had invariant masses between 2.8 
and 3.4 GeV/c 2 . The Level-2 J/tp trigger accepted 19 
of these events, yielding an estimate of 0.73±o;°o (stat), 
in very good agreement with the trigger simulation re- 
sult. These results verify that the trigger performance is 
similar for real data and simulations. 



The J/tp acceptance and efficiency e a 



cc—eff 



is deter- 



mined with a GEANT based Monte Carlo simulation 
of the PHENIX experiment. The detector response has 
been tuned to reproduce the resolution and performance 
of the real detector. The efficiency includes not only the 
tracking efficiency, but also the probability for passing all 
of the electron identification selection cuts. The electron 
identification efficiency determined by the Monte Carlo is 
verified by a clean electron sample from conversion pho- 
tons. We also account for run by run efficiency changes by 
counting the relative number of reconstructed electrons 
and positrons per event in our data sample. We show 
the PHENIX acceptance and efficiency as a function of 
transverse momentum in Figure |3J 

Since we do not have the statistics to determine the 
transverse momentum distribution of the J/tp, we must 
employ a model for the px dependence to determine an 
overall acceptance and efficiency. We use two different 
functional forms for the px distributions to test the model 
sensitivity of our acceptance. We use an exponential in 
Pt and an exponential in p T as motivated by fits to J/tp 
data at lower energies [22 ■ The two models give similar 
acceptance values given a common < px > value input. 
The largest uncertainty comes from the value of < px > 
assumed. PHENIX has measured J/tp production in 
proton-proton reactions at t/s — 200 GeV and finds a 
<PT > = l.80±0.23 (stat) ±0.16 (sys) GeV/c Q. We 
use this value to determine our acceptance and efficiency 
averaged over all px- The J/ip < px > in Au-Au col- 
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FIG. 3: The PHENIX J/tp acceptance x efficiency as a func- 
tion of the J/tp transverse momentum is shown. Most of the 
acceptance is with one lepton into each of the two PHENIX 
central spectrometers. This contribution peaks at pr — 
and decreases with increasing pr- The rise in the accep- 
tance at high pr is from contributions where both electron 
and positron are accepted into one of the PHENIX central 
spectrometers. 



TABLE H: The J/tp acceptance x efficiency and the central- 
ity dependent efficiency are shown for three exclusive Au-Au 
centrality event classes. 



Centrality bin 


^QrCC — £ f f 


tcent 


00-20% 


0.0027 ±K (sys) 


0.61 ±0.06 (sys) 


20-40% 




0.78 ±0.08 (sys) 


40-90% 


ft 


0.90 ±0.09 (sys) 



lisions may differ from that in proton-proton reactions. 
Therefore we vary the < px > from 1.0 to 3.0 GeV/c to 
determine our model dependent systematic errors. We 
assume that the J/tp rapidity distribution is flat over the 
range —0.35 < y < 0.35 where we measure. The final 
value for the J/tp acceptance and efficiency is shown in 
Table ILT1 This acceptance and efficiency has a 20% sys- 
tematic error from uncertainties in matching the Monte 
Carlo to the detector response, a 10% systematic error 
from run-to-run variation corrections, and a ±24% sys- 
tematic error from the uncertainty in the < px >■ 

Our tracking and electron identification efficiencies ex- 
hibit a centrality dependence due to overlapping hits and 
energy contamination in the calorimeter. We determine 
this dependence by embedding Monte Carlo J / tp into real 
data events of different centrality selections. The corre- 
sponding efficiency factor e cent varies from 56% for the 
0-5% most central events to 98% for the 85-90% most 
peripheral events. 

The final values for the embedding efficiency in our 
wide centrality bins are sensitive to the true centrality de- 
pendence of the J/tp production. In order to estimate the 
systematic error due to this uncertainty we assume two 
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different centrality dependence models: (1) binary colli- 
sion scaling and (2) participant collision scaling. Within 
our centrality ranges, we find that these two models yield 
less than a 5% difference and we include this in our sys- 
tematic error. We assign an additional 10% systematic 
error to account for uncertainties in the Monte Carlo em- 
bedding procedure. The centrality dependent efficiency 
values are shown in Table ITU 

In our B dN/dy calculation, we have added the system- 
atic errors from all of the contributing factors in quadra- 
ture and find +35% and -41% total systematic error on 
the invariant yield in each of the centrality ranges. The 
dominant systematic error results from the uncertainty 
in the < pr > of the J/ip distribution. 

VI. RESULTS 

The B dN / dy\ y =a values for the three exclusive cen- 
trality selections are shown in Table IIIII We have cal- 
culated using a Glauber model [2l| the number of ex- 
pected participating nucleons N par t and the number of 



expected binary collisions N co u for each centrality range. 
These results are shown in Table Hvl in addition to the 
B dN/dy\ y=0 values divided by the expected number of 
binary collisions. 



TABLE III: We show the statistically most likely J/ip invari- 
ant yield (B dN/dy\ y =o) value and the 68% confidence in- 
terval for peripheral (40-90%) and mid-central (20-40%) col- 
lisions. We also show the 90% confidence level upper limit 
and the systematic error on this limit for all three different 
centrality ranges of Au-Au collisions. 





B dN/dy\ v=0 


(xlO" 4 ) 


Centrality 


Most Likely Value 


90% C.L.U.L. 


00-20% 


N.A. 


6.08+1.56 (sys) 


20-40% 
40-90% 


4.001^ (stat)+_{f (sys) 
0.86lS;S (stat)±° f 5 (sys) 


7.19 + 2.43 (sys) 
1.60 + 0.54 (sys) 



The PHENIX result for the J/ip invariant yield in 
proton-proton induced reactions at ^/s=200 GeV at mid- 
rapidity [24| is 



B dN/dy\ y=0 (pp) = 1.46 ± 0.23 (stat) ± 0.22 (sys) ± 0.15 (abs) x 10~ 6 (4) 

I 



The systematic error (abs) represents the uncertainty of 
the normalization of the total proton-proton invariant 
yield. 

We show in Figure 01 the results from the three Au-Au 
centrality bins and the proton-proton data normalized 
per binary collision as a function of the number of partic- 
ipating nucleons. Note that for proton-proton reactions, 
there are two participating nucleons and one binary col- 
lision. 



VII. DISCUSSION 

Despite the limited statistical significance and system- 
atic uncertainty of these first J/ip results, we can address 
some important physics questions raised by the numerous 
theoretical frameworks in which J/ip rates are calculated. 

We show in Figure binary scaling expectations as 
a gray band. We also show a calculation of the sup- 
pression expected from "normal" nuclear absorption us- 
ing a o C c-N = 4.4 mb H3 and 7.1 mb @, El- A re- 
cent measurement in proton-nucleus collisions at lower 
energies favors the smaller absorption cross sec- 
tion, thus underscoring the importance of measuring J/ip 
in proton(deuteron)-nucleus collisions at RHIC energies. 
We also show the NA50 suppression pattern relative to 
binary scaling |3j , normalized to match our proton-proton 
data point at 200 GeV. The data disfavor binary scaling, 



while they are consistent with "normal" nuclear absorp- 
tion alone and also the NA50 suppression pattern mea- 
sured at lower energies, within our large statistical errors. 

One model calculation ^(| including just the "normal" 
nuclear and plasma absorption components at RHIC en- 
ergies is shown in Figure The higher temperature (T) 
and longer time duration of the system at RHIC lead to 
a predicted larger suppression of J/ip relative to binary 
collision scaling. This s pec ific model and in general 
this class of models |§l l32j. cannot be ruled out at this 
time due to our null result (90% confidence level upper 
limit) for the most central collisions. 

Many recent theoretical calculations also include the 
possibility for additional late stage re-creation or coales- 
cence of J/ip states. In [l||, they include both break-up 
and creation reactions D + D<-+J/ip + X. At the lower 
fixed target CERN energies, this represents a very small 
contribution due to the small charm production cross sec- 
tion. However, at RHIC energies, where in central Au-Au 
collisions there are of order 10 cc pairs produced, the con- 
tribution is significant. The sum of the initial production, 
absorption, and re-creation as shown in Figure is also 
consistent with our experimental data. 

A different calculation ^{| assumes the formation of 
a quark-gluon plasma in which the mobility of heavy 
quarks in the deconfined region leads to increased cc co- 
alescence. This leads to a very large enhancement of 
J/ip production at RHIC energies for the most central 
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TABLE IV: We show the number of participating nucleons and the number of binary collisions for three different centrality 
ranges of Au-Au collisions, and the associated systematic errors. We show the statistically most likely value for the J/ip 
invariant yield (B dN/dy\ y =a) divided by the expected number of binary collisions for peripheral (40-90%) and mid-central 
(20-40%) collisions. We also show the 90% confidence level upper limit and the systematic error on this limit for all three 
different centrality ranges of Au-Au collisions. The systematic error in the invariant yield per binary collision does not include 
the systematic error in the expected number of binary collisions. This error contribution is negligible for the central and 
mid-central categories and would increase the systematic error for the peripheral category by 6%. 



Centrality 


Npart 


Nooll 


B dN / dy\ y= o per 
Most Likely Value 


binary collision (xlO 6 ) 

90% C.L.U.L. 


00-20% 


280 ±4 


779 ± 75 


N.A. 


0.78+0.20 (sys) 


20-40% 
40-90% 


140 ±5 
34 ±3 


296 ± 31 
45 ±7 


1.351S;™ (stat)±° :lt (*V« 
1.9llS:5? (stat)±° ( £ (sys 


) 2.43 + 0.82 (sys) 
) 3.55+1.21 (sys) 
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FIG. 4: The J/ip invariant yield per binary collision is 
shown for proton-proton reactions and three exclusive central- 
ity ranges of Au-Au reactions all at ^/snn = 200 GeV. For 
the proton-proton reactions, we show the most likely value 
as a data point (square), the statistical error, and the esti- 
mated systematic errors as brackets. For the three Au-Au 
data points, we show as arrows the 90% confidence level up- 
per limits. The bracket above the limit includes the estimated 
systematic error on these limits. In the case of the periph- 
eral and mid-central ranges, we also show, as a square marker, 
the statistically most likely value and as two horizontal dashes 
the 68% confidence interval. The gray band indicates binary 
scaling and the width is the quadrature sum of the statistical 
and systematic error on our proton-proton data point. For 
the Au-Au points, the systematic error in the invariant yield 
per binary collision does not include the systematic error in 
the expected number of binary collisions. This error contri- 
bution is negligible for the central and mid-central categories 
and would increase the systematic error for the peripheral 
category by 6%. 



reactions. The model considers the plasma temperature 
(T) and the rapidity width (Ay) of charm quark pro- 
duction as input parameters. Shown in Figure [H] are the 
calculation results for T = 400 MeV and Ay = 1.0, 2.0, 
3.0, 4.0. The narrower the rapidity window in which all 
charm quarks reside, the larger the probability for J/ip 
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FIG. 5: (Color online) The J/ip invariant yield per binary 
collision is shown from proton-proton reactions and three ex- 
clusive centrality ranges of Au-Au reactions all at ^/sjvjv = 
200 GeV. The lines are the theoretical expectations from "nor- 
mal" nuclear absorption with a C c-N = 4.4 mb (solid curve) 
and 7.1 mb (dashed curve) cross section. The stars are the 
J/ip per binary collision measured by the NA50 experiment 
at lower collision energy. In order to compare the shapes of 
the distribution, we have normalized the NA50 data to match 
the central value for our proton-proton results. 



formation. Ay — 1.0 is consistent with the three dimen- 
sional spherically symmetric thermal distribution, and 
results in a charm yield at midrapidity that is inconsis- 
tent with the PHENIX preliminary charm yield as deter- 
mined from single electron measurements [l4|. Ay = 4.0 
is consistent with expectations from factorized QCD and 
PYTHIA with CTEQ5L structure functions [li|. All of 
these parameters within this model predict a J/tp en- 
hancement relative to binary collisions scaling, which is 
disfavored by our data. 

Another framework for determining quarkonia yields 
is to assume a statistical distribution of charm quarks 
that may then form quarkonia. A calculation assuming 
thermal, but not chemical equilibration |17| is shown in 
Figure and is also consistent with our data. 
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FIG. 6: (Color online) The J/ip invariant yield per binary 
collision is shown from proton-proton reactions and three ex- 
clusive centrality ranges of Au-Au reactions all at ^snn = 
200 GeV. The lowest curve is a calculation including "normal" 
nuclear absorption in addition to substantial absorption in a 
high temperature quark-gluon plasma |lol . The curve above 
this is including backward reactions that recreate J/ip. The 
statistical model (l7j result is shown as a dotted curve for mid- 
central to central collisions just above that. The four highest 
dashed curves are from the plasma coalescence model |l5f| for 
a temperature parameter of T = 400 MeV and charm rapid- 
ity widths of Ay = 1.0, 2.0, 3.0, 4.0, from the highest to the 
lowest curve respectively. 

Significantly larger data sets are required to address 
the various models that are still consistent with our first 
measurement. Key tests will be the pr and xf depen- 
dence of the J/ip> yields, and how these compare with 
other quarkonium states such as the xp'. 

VIII. SUMMARY 

PHENIX has shown first results on J/ ip production in 
Au-Au collisions at ^/snn = 200 GeV at mid-rapidity 
as measured via electron-positron pairs. We find that 
models that predict J/tp enhancement relative to binary 



collision scaling are disfavored, while we cannot discrim- 
inate between various scenarios leading to suppression 
relative to binary scaling. 

This first measurement from PHENIX will be followed 
with high statistics measurements in both the electron 
channel at midrapidity and at forward and backward ra- 
pidities in the PHENIX muon spectrometers. Such mea- 
surements are expected in the next few years and will 
address the full range of heavy quarkonia production and 
evolution models. 
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